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Abstract 

We study zero temperature phase transitions in two classes of random quan- 
tum systems -the g-state quantum Potts and clock models. For models with 
purely ferromagnetic interactions in one dimension, we show that for strong 
randomness there is a second order transition with critical properties that 
can be determined exactly by use of an RG procedure. Somewhat surpris- 
ingly, the critical behaviour is completely independent of q (for 2 < q < oo). 
For the q > 4 clock model, we suggest the existence of a novel multicritical 
point at intermediate randomness. We also consider the T = transition 
from a paramagnet to a spin glass in an infinite range model. Assuming that 
the transition is second order, we solve for the critical behaviour and find q 
independent exponents. 
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The effects of randomness on the properties of quantum many-body systems have been 
the subject of experimental and theoretical studies for many years [|1|. Despite this, un- 
derstanding of phenomena where quantum effects, randomness, and interactions are all 
important is rather poor, and there are as yet few reliable theoretical techniques. Recent 
work has focussed attention on simple quantum statistical mechanical systems with 
randomness as a useful starting point to obtain insights into such phenomena, and some 
progress has been made. In this paper we study analytically the effects of disorder on the 
properties of two classes of quantum models with discrete symmetry. These are to be re- 
garded as quantum versions of the classical g-state Potts and the g-state clock models [§. 
Studies of the corresponding classical models have yielded a fair amount of insight into the 
competing effects of randomness, interactions, and thermal fluctuations [0. For instance, 
classical Potts spin glasses have been studied extensively as a paradigm for understanding 
the properties of orientational glasses . 

We first consider the zero temperature quantum phase transition in random (purely fer- 
romagnetic) 1-dimensional quantum Potts and clock chains. For the 1-d random transverse 
field Ising model (RTFIM), which (as we shall see below) is the q = 2 case of these models, 
a wealth of essentially exact information has recently been obtained in a remarkable paper 
by Fisher using a real space renormalization group procedure. Here we show that this 
procedure can similarly be used to obtain exact critical properties of strongly random g-state 
Potts and clock chains for all g, and implies, remarkably, that there is no g dependance in 
any of the exponents or the scaling functions. This is in stark contrast to the pure problem 
where the properties of the transition depend crucially on the value of g In addition, 

considerations on the effects of weak randomness suggest the possibility that for any amount 
of randomness, the Potts model for any g and the clock model for 2 < g < 4 are described 
by the strong randomness fixed point. For the clock chain with g > 4, we suggest the exis- 
tence of a multicritical point at a finite strength of randomness. Next we consider the zero 
temperature transition from a paramagnet to a spin glass in infinite-ranged quantum Potts 
and clock models. Again building on work on the g = 2 case |^,^, and assuming a second 
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order transition, we find that the critical properties are independent of q. 

The models are defined in terms of a variable that can assume q possible states (which 
we denote |0) , |1) , ■ ■ ■ , |g — 1)) on the sites of a d-dimensional lattice. The classical Potts 
(clock) interaction in the presence of a uniform external "magnetic" field H along the '0' 
direction is 

np,int = - E Jv^n^n, " 2if ^('^".,0 ' ') {Potts) (l) 

{i,3) i ^ 

T-icint = - E '^■^ij — ('^i -"^j)) - '^H^cos (Clock) (2) 

{ij) ^ i ^ 

We introduce quantum fluctuations into these models by adding at each site a "transverse 
field" term that attempts to change the state of the variable at that site. Thus we consider 
the quantum Hamiltonians 

np = -J2hj E -l^i) + ^P,int (Potts) (3) 
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= - E H\n^){n. + 1| + h.c) + nc,^nt (Clock) (4) 

\n,=0 / 

(We identify \ni + q) = \ni)). Through out the paper we will assume that the hi and Jij are 
independent random variables drawn from some distributions Pi(h) and P2(J) respectively. 
Note that at = 0, the Hamiltonian Tip is invariant under a global permutation \n) \n') 
of the states at each site. For Tic, the symmetry is a global cyclic rotation \n) —>■ \n + 1). 
Clearly for q = 2, both these models reduce to the transverse field Ising model. For general 
q, just as in the Ising case, the "transverse field" term plays the role of a kinetic energy 
that opposes the tendency to order due to the interaction term. Also as in the Ising case, 
the d-dimensional g-state quantum Potts (clock) model Eqn. ^ (^) at zero temperature may 
be regarded as the transfer matrix in the r-continuum limit of a d + 1-dimensional g-state 
classical Potts (clock) model with disorder constant along one direction. 

In the absence of disorder, the mapping to the classical d + 1-dimensional pure problem 
provides a rather complete picture of the possible phases and the transitions between them. 
For instance (at zero H), the ferromagnetic (i.e J > 0, h > 0) quantum Potts chain has 



a first order transition for q > 4, and a second order transition for g < 4 (for which all 
the exponents are known exactly and depend on the value of q) P,p^. The ferromagnetic 



clock chains, on the other hand, have, for g > 4, a quasi-long-range ordered (QLRO) phase 
sandwiched between a truely long-range ordered phase and a disordered phase [§] . For 
g < 4, the quasi- long-range ordered phase disappears and is replaced by an ordinary second 
order phase transition for which again all the exponents are known exactly Below we 

will see that randomness drastically modifies this picture. 

We start with 1-d chains in which all the /ij's and Ji/s are random but positive and 
H = 0. Defining the total magnetization as M = with Si = 6^,0 — ^ for the Potts 

model and Si = cos for the clock model, it is clear that as the overall relative strength 
of the /ij's is decreased, there will be a transition from a phase with M = to one with 
M 7^ 0. We assume that the randomness is strong and follow closely a real-space RG 
procedure used by Fisher p[ to extract an enormous amount of information on the RTFIM 



(g = 2 case). The basic idea behind this procedure |Tl]] is to successively eliminate the 
strongest coupling Q = max{hi, Jj,i+i} in the chain and get an effective Hamiltonian for the 
low energy degrees of freedom. First consider the case when the maximum coupling is a 
field, say h^. We eliminate the site i, and obtain, using second order perturbation theory, 
a new effective bond between the sites i — 1 and i + 1 of strength J = (Ji-i^Ji^i+i) / nhi 
where k is g/2 for the Potts model and (1 — cos(^))/(l + 6q^2) for the clock model. On 
the other hand, if the maximum coupling is a bond we replace the sites i and i + 1 

by a single Potts (or clock) degree of freedom with an effective field h = (hihi^i) / nJi^i+i 
where k is the same as before (as may be expected from a duality which these models can be 



shown to possess [[l^]). The g dependance is only through n. As in Ref. 0, we convert these 
recursion relations into flow equations for the distributions Pi(C = In(i7/J),r = ln{Qj/Q)) 
and P2{l3 = ln{Q/h),r = ln{Qj/Q)) where Qj is the initial value of the maximum coupling 
and find 

dPi dPi 



dT dC 



Pi(Pi(o,r)-P2(o,r)) + 



+ P2(o, r) J dCidC2Pi{Cu r)Pi(C2, r)5(c - Ci - C2 - in^) (5) 

and similarly with Pi ^ P2 (as expected from duality). When now we rescale and look 
for critical fixed points, the value of In k becomes irrelevant at low energies (so long as it 
is finite). The resulting probability distributions in the scaling limit described by the fixed 
point are independent of q. In general it is necessary to keep track of two joint distributions 
- that of bond lengths and bond strengths at scale F and that of cluster lengths, their 
magnetic moments and "transverse" field strengths at scale F. Both these distributions will 
be independent of q in the scaling limit. Thus the value of q merely determines a high-energy 
cutoff Ec ~ i^/e^''"^"' (below which one should be in order to observe scaling behaviour) but 
does not affect the scaling behaviour itself. Note that since this cutoff goes to zero as g — > 00, 
our results hold only for finite q. 

At this point it is tempting to conclude that all the Potts and clock chains will have 
identical critical properties. However we note that, in general, identical probability distribu- 
tions do not necessarily imply identical physical properties. For instance, the pure problems 
trivially have identical distributions but have very different properties. However, as shown 
by Fisher |^ , the distributions of the logarithmic couplings ( and (3 become infinitely broad 
asymptotically at low energies. It is then straightforward to see that due to this extreme 
broadness all physical quantities (such as magnetization or mean correlation function) are 
described by g-independent scaling functions and exponents. The only g-dependence occurs 
in some non-universal constants. We illustrate this point with the example of the scaling 
of the magnetization M{H, 6) as a function of small external applied field H{> 0) and a 



dimensionless measure 5 of the deviation from criticality [13]. In the presence of a magnetic 
field H, the energy levels of an otherwise-free cluster of magnetic moment fi split into a 
ground state \n = 0) and other excited states with a gap Eh = 2fiqH for the Potts case and 
(atleast) Eh = 2jj,H{l — cos ^) for the clock case. Proceeding exactly as for the RTFIM, 
we stop the RG when the maximum coupling ^7 ~ Eh- Due to the extreme broadness of 
the distribution, an asymptotically exact expression for M{H, S) is obtained by aligning all 
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the remaining clusters in the direction of the magnetic field. Thus 

M {H , 5) = ft X (total number of active (i.e undecimated) spins at scale F/^ = In — — 

H 

+ corrections 

where /i and Dh are non-universal and possibly q dependent constants. The key point 
now is that the number of active spins at a given scale F is entirely a property of the 
joint distribution of cluster lengths, magnetic moments, and field strengths which is q- 
independent in the scaling limit. Consequently, the universal scaling function describing 
M{H, 6) will be the same for all q. The q dependence is only in the non-universal quantities 
/i and Dh- Through similar reasoning, one can establish that the mean correlation function 
C{x) = (sqSx) is also described by a g-independent scaling function. 

Thus the Potts and clock chains for any q (with strong ferromagnetic randomness) do 
indeed have critical properties identical to those of the RTFIM (the q = 2 case) for which 
detailed results are available . We point out some salient features below. The spontaneous 
magnetization vanishes at the transition with exponent /? = ■ The mean and typical 
correlation functions at the critial point decay as C{x) ~ and — In Ctyp{x) ~ \f\x\ re- 
spectively for large In the disordered side, there are 2 correlation lengths, characterizing 
the decay of mean and typical correlations, which diverge with exponents z/ = 2 and 9 = 1 
respectively. The magnetization scaling function is known exactly, and the scaling function 
for the mean correlation function known up to the solution of a linear ordinary differential 
equation. Asymptotically at the critical point "lengths" scale as the square of the logarithm 
of "energies" (unlike most other quantum transitions where lengths scale as a power of the 
energies). On either side of the transition there are Griffiths regions. In the Griffiths region 
of the disordered phase, the order-parameter susceptibility diverges as T ^ as a power 
with an exponent weaker than Curie. Throughout this region, the magnetization increases 
as a power (with logarithmic corrections) of an applied external magnetic field with a con- 
tinuously varying exponent. Similarly in the Griffiths region of the ordered side the stiffness 
vanishes for an infinite system, and the susceptibility diverges as T — as a power with an 
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exponent that is stronger than Curie. Very far from the critical point, there are of course 
the more conventional strongly ordered and disordered phases. 

The RG procedure used above is valid only if the randomness in the initial distributions 
Pi and P2 is strong. Clearly it cannot address the question whether, if the initial distribution 
is narrow, the low energy properties will still be described by the strong randomness fixed 
point found above. Some insights on this matter are provided by considering the effect 
of weak randomness on the pure systems. So long as the pure transition is second order, 
weak randomness is relevant at the fixed point ii u < 2/d (the generalized Harris criterion 
|1^) where d is the spatial dimension. From the known values of u |]6|JT0|] for g < 4 for the 



Potts and clock chains, we conclude that weak randomness is indeed relevant. The simplest 
scenario then is that for g < 4, the RG flows take the system to the strong randomness fixed 
point for any amount of randomness in the initial distributions. The situation however is 
different for g > 4. We discuss the Potts and clock cases separately. The pure Potts chain 
for g > 4 has a first order transition. For first order phase transitions in classical systems, it 
has been argued that any amount of "bond" randomness converts the transition to second 
order |TH| for d < 2. Extension of this argument to quantum systems would then suggest 



that the random g > 4 quantum Potts chain has a second order transition for any amount 
of randomness. We conjecture that this transition is described by the strong randomness 
fixed point. 

We now turn to the clock model for g > 4 where there is an intermediate quasi-long-range 
ordered phase separating a truely ordered phase from the paramagnetic phase. The phase 
transitions into the QLRO phase from either side are of the Kosterlitz-Thouless type with 
u = 00, therefore implying by the Harris criterion that weak randomness is irrelevant. The 
QLRO phase is described by a line of fixed points. A straightforward perturbative calculation 



shows that weak randomness is irrelevant along this entire line |T^ . Thus weak disorder does 
not change the nature of the phase diagram in this case while for strong disorder, as we have 
seen earlier, there is a single second order phase transition (surrounded by Griffiths regions) 
and no QLRO phase. Understanding how the phase diagram changes as the strength of 



the disorder increases is an interesting open question. We speculate that as the strength 
of the disorder is increased, the two hues of Kosterhtz-Thouless transitions bounding the 
QLRO phase merge at a muhicritical point. Beyond this point the QLRO phase disappears, 
and there is a single second order transition described by the strong randomness fixed point 
(See Figure |I]). We do not however have any strong arguments ruling out more complicated 
scenarios (such as for instance, a new intermediate phase separating the QLRO phase from 
the region with a single second order transition). 

Having described the 1-dimensional random ferromagnetic systems in some detail, it 
is natural to ask if there are other non-trivial yet solvable cases. One such example is 
provided by infinite-range spin glass models. As is well-known, classical SG models display 
a complicated and rich structure even for infinite-range interactions. There has been some 
recent progress in understanding the zero temperature quantum phase transition into a spin 
glass phase in the transverse field Ising model (the q = 2 case of our models) . The basic 
idea is quite straightforward |[T^. One first performs the disorder average using the replica 
trick and reduces the problem to an effective single site problem with a self-consistency 
condition on the auto-correlation function. At the critical point it is possible to solve the 
self-consistency condition by making a suitable ansatz for the long-time behaviour of the 
auto-correlation function. For the Ising case, one finds that the critical auto-correlation 
decays as at large imaginary time r P,^. In addition, in the paramagnetic side, there 
is an energy gap that vanishes on approaching the transition with an exponent 1/2 with 
logarithmic corrections. Repeating the analysis for the g-state quantum Potts or clock spin 
glass model, and assuming a second order transition, we find [0] that the critical exponents 



once again remain the same as the Ising case. However the assumption of a second order 
transition may be questionable (atleast for the Potts and odd-g clock models) since the 
finite temperature phase transition from the paramagnetic side in the corresponding classical 



models is not of the conventional second order type |T8 . 

In summary, we have shown that for strongly random ferromagnetic quantum g-state 
Potts and clock chains, the critical properties are completely independent of q for 2 < q < oo. 
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This result is a priori surprising as in the pure models the possible phases and the transitions 
between them are known to depend crucially on the value of q. We have also studied 
the zero temperature transition from a paramagnet to a spin glass in these models with 
infinite range interactions. If the transition is second order, the critical exponents have 
g-independent values. We conclude by noting the following implications of our work and by 
raising some open questions. Our results imply that for classical 2-d Potts and clock models 
with disorder correlated along one direction, the critical properties of the finite temperature 
phase transition are independent of the value of q. It is interesting that numerical simulations 
by Chen et. al. |jl9l on the classical 2-d 8-state Potts model with uncorrelated disorder find a 
second order transition with exponents equal to the classical 2-d Ising values. This suggests 
the possibility that the g-independence found here with correlated disorder also holds with 
uncorrelated disorder. However an analytic calculation by expanding in g — 2 |2^ does 
find g-dependence in the exponents. Further studies to address this issue will be welcome. 
From our discussion of the phase diagram for the clock chains with g > 4 (Figure 0), there 
arises the possibility of a novel multicritical point at finite randomness. Verification of the 
existence of this point, perhaps numerically, is an interesting open problem. Also interesting 
is the question of whether this g independence persists in higher dimensional ferromagnetic 
models. For the spin glass models, one approach to go beyond our results, atleast in high 
enough dimensions, would be to study a Landau theory analogous to that developed for the 
Ising and rotor models . Such a study may also resolve the issue of whether the transition 
is second order or not. 

We thank Subir Sachdev, R.Shankar, N.Read, and D.S.Fisher for useful discussions and 
comments. This research was supported by NSF Grants No. DMR-92-24290 and DMR-91- 
20525. 
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FIG. 1. One possible (schematic) phase diagram for the q > 4 random quantum clock chain. 
g = \nh is a parameter that measures the strength of quantum fluctuations and V is & measure of 
the strength of the randomness. 
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